The CMB bispectrum from recombination 
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We compute the cosmic microwave background temperature bispectrum generated by nonlin- 
earities at recombination on all scales. We use CosmoLib2"'^, a numerical Boltzmann code at 
second-order to compute CMB bispectra on the full sky. We consistently include all effects except 
gravitational lensing, which can be added to our result using standard methods. The bispectrum 
is peaked on squeezed triangles and agrees with the analytic approximation in the squeezed limit 
at the few per cent level for all the scales where this is applicable. On smaller scales, we recover 
previous results on perturbed recombination. For cosmic-variance limited data to /max = 2000, its 
signal-to-noise is S/N = 0.47 and will bias a local signal by /^l — 0.82. 



Introduction: The non-Gaussianity of cosmological 
perturbations is a powerful discriminator of early uni- 
verse models. In particular, a future detection of local- 
type non-Gaussianity in the forthcoming Planck data [T] 
will rule out all single-field models with an adiabatic 
attractor and Bunch-Davies vacuum. However, nonlin- 
earities connecting the initial conditions to the observed 
CMB can potentially contaminate the primordial signal 
and must be fully understood. The most important of 
these effects are at late-time, such as integrated Sachs- 
Wolfe (ISW)-lensing correlation [2l|3|, but second-order 
sources during recombination may also be relevant. 

Computing the full bispectrum from recombination re- 
quires to resort to a numerical Boltzmann code. Al- 
though there has been an intense effort to derive the 
complete second-order equations [IHS], getting the bis- 
pectrum is very complicated and complete control of all 
effects has not been reached yet. Currently, the only 
complete numerical calculation is from the code CM- 
Bquick [2, which claims a contamination of /^l ^ 5 
from recombination j7|. However, this is in tension with 
partial calculations focussing on particular limits. For 
super-horizon modes at recombination the contamina- 
tion is comparable to /fl^L ^ —1/6 (TUl HJ- On sub- 
Hubble scales two effects have been computed: fluctu- 
ations in the free-electron density, inducing local-type 
non-Gaussianity by delaying the time of recombination 
[121 113| , and the nonlinear evolution of dark matter in- 
ducing equilateral-type non-Gaussianity in the photon 
temperature fluctuations [HI [T5] . None of these effects 
will relevantly bias a local signal detection. 

Given the current status and the relevance of the re- 
sult, it is important to shed some clarity on the full bis- 
pectrum from recombination, estimate its contamination 
and detectability. In this letter we compute the full-sky 
bispectrum from recombination on all scales. We use 
CosmoLib2"'^, a second-order numerical Boltzmann code 
which will be made publicly available in the future. 

Fortunately, we have a robust theoretical guidance 
which allows to put the code on solid grounds. The bis- 
pectrum generated at recombination can be computed 



analytically in the squeezed limit and when one of the 
three modes is outside the horizon at recombination [18J- 
[2T]. The argument is analogous to that of the single- 
field consistency relation [52]. After inflation, the effect 
of a super-horizon mode can be traded with a scale re- 
definition on the short ones. Thus, for scale-invariant 
perturbations three modes of the primordial curvature 
perturbation C, are decorrelated in the squeezed limit. Of 
course, this is true only if C, appears exponentially in 
the spatial metric. Analogously, since T (x 1/a the bis- 
pectrum of three super-horizon modes of O defined by 
T — Te® vanishes in the squeezed limit. The nonlin- 
ear relation converting the "Gaussian" variable 8 to the 
usual temperature fluctuation O = ST/T gives a local 
modulation inducing a temperature reduced bispectrum 

hihh = + Ci^Ci^ + Ci^Ci^ -f fcii/aia , (1) 

where Ci is the temperature spectrum and bi-^i^i^ is the 
reduced bispectrum of O. Moreover, local physics during 
recombination is unaffected by the presence of a super- 
horizon mode of C, which just acts as an unobservable 
scale redefinition. However, when the long mode eventu- 
ally re-enters the horizon, it induces an angular rescaling 
of the observed 2-point function in different regions in the 
sky. As the temperature spectrum substantially departs 
from scale-invariance due to oscillations in the photon- 
baryon plasma, this induces a temperature bispectrum 
for O given by 

cr 

where Cj is the cross-correlation spectrum between O 
and C at recombination and li <C ^2.3- The symmetriza- 
tion guarantees that corrections to this expression are of 
order C'(/i/?2,3)^- Corrections due to sub-Hubble physics 
scale as oIIi/IhY, with Ih = (iJa/c^ )rcc?7o - HO [II]. 

The local redefinition ([T]) and the spacial rescaling ^ 
are not the only contributions to the bispectrum in the 
squeezed limit. Lensing close to the last scattering sur- 
face displaces the 2-point function inducing an isotropic 
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and anisotropic rescaling of tfie coordinates [31 [11] [15] . 
Altfiougli in our numerical calculation we do not include 
lensing, the separation is clear: eqs. ([T]) and ^ are the 
analytic approximation in the squeezed limit to the bis- 
pectrum computed by integrating the photon distribu- 
tion along unperturbed geodesies. 

The code: CosmoLib2"'^ is a Boltzmann code based 
on CosmoLib [T7]. It solves Einstein and Boltzmann 
equations in Poisson gauge for all cosmological species 
up to second-order in the perturbations and computes the 
full-sky CMB bispectrum by integrating the second-order 
photon brightness along the line of sight and projecting 
it into harmonic space. 

For the Boltzmann solver we have employed the equa- 
tions given in Ref. [7 complementing them with those 
derived in several other references [IHSl [E]j paying at- 
tention to correct typos. The initial conditions are set 
on super-horizon scales assuming a minimal model of in- 
flation where the primordial curvature perturbation ^ is 
approximately Gaussian. More generic adiabatic initial 
conditions are possible. Perturbed recombination is in- 
cluded using DRECFAST [23]. The perturbation in the 
free-electron density perfectly agrees with that computed 
by [6] on large and small scales; < 10% disagreement is 
registered on scales of order Hubble at recombination, 
due to the inconsistent treatment of metric perturbations 
in DRECFAST. This wiU be corrected in the future. 

Second-order solutions with one of the modes super- 
horizon computed by the code agree with those obtained 
by a coordinate transformation as explained in |19| . On 
small scales and matter dominance, they agree with the 
analytic solutions that are known in the literature |24j . 
Moreover, we have extensively used CMBquick to com- 
pare our Boltzmann solutions. Improvements with re- 
spect to such a code are: more precise treatment of 
the tight-coupling approximation necessary to evolve the 
early-time solutions from the initial conditions; different 
line-of-sight integration which reduces the errors due to 
truncations in the niultipole expansion of the source, as 
described below; bispectrum computed for the full sky 
and with a different integration method. 

Bispectrum computation: The photon brightness 
is defined from the photon one-particle distribution / as 
I{r],x,n) = J dpp'^ fiji^x^p), where is the momentum 
of photons in the local inertial frame and = pn*, riin^ = 
1. We use the metric in Poisson gauge, goo = — a^(?7)e^*, 
goi = a^{v)^i and = a?{vi){e~'^'^ 5ij +Xi])i with uji^i = 
and Xii = = Xij,j- From the Boltzmann equation 
one can write an equation for the fractional brightness 
A ^ 51/1 



d 
drj 



(A + 4$)-4A(*-$,,:n')-£; = ~t{1 + 5^ + ^)F , (3) 



where d/drj is the advective derivative, a dot denotes 
the derivative with respect to rj and f is the unper- 
turbed differential optical depth. On the l.h.s., E = 



4(i' -I- — Aujin^ — 2xijn^n^ is the redshift in photon's 
energy due to integrated effects: the ISW contribution 
(whose second-order part is the Rees-Sciama (RS) effect), 
and the vector and tensor contributions, respectively [llj . 
The r.h.s. comes from the collision term of the Boltz- 
mann equation: F can be read off, for instance, from the 
r.h.s. of eq. (78) of [5], with the notation for $ and ^ 
interchanged; Se is the perturbation of the free-electron 
density. 

The second term on the l.h.s. of (|3| is problematic: it 
must be integrated along the line of sight and couples 
all multipole moments of the fractional brightness to the 
gravitational potentials. A crucial simplification comes 
from trading this term for a boundary term and one pro- 
portional to T, contributing only at recombination. By 
replacing <I>,in' by d^/drj — <i> and dividing eq. ([s]) by 
1 -|- A we can rewrite it, up to second order, as 



^-/™ d>,'f(l+5.+<t.)_ 



d \/~ 



di] 



(a + 4$) e-^^'^'^'^d 



+<5e+*) 



(4) 



-E = -tR = -f(l + 6e + ^){F + A + 4$) 



We have introduced the variable A = A — A^/2, defined 
from the brightness in exponential form, / ~ le^, and 
F = F(l — A). As explained above, A is Gaussian in the 
squeezed limit and on super-horizon scales. Analogously, 
one can show that F on the r.h.s. of Q becomes Gaussian 
in the same limit \W\ . Integrating this equation formally 
yields 

rvo 

A(77o, f = 0, 7i) = / d77e~^('')5(7?, f (r;), n) , (5) 
Jo 

S={E-tR)(^1+ df]'T {6e + $) j , (6) 

where S is computed along photon trajectories. Here we 
integrate along a straight trajectory x{r]) — fi{ri — rjQ). In 
doing this we are neglecting the Shapiro time-delay and 
gravitational lensing effects. While the former can be 
safely neglected [25], gravitational lensing can be taken 
into account using standard treatments [21 [5] . Consistent 
inclusion of lensing and time-delay effects in the code is 
left for future work. 

The observed brightness A (770, fi) in ([s]) can be decom- 
posed into spherical harmonics. 



6(770, n) = ^A(7;o,7^) = ^di„-Xim{n) , 

Im 



(7) 



where the factor of 1/4 in the first equality comes from 
the relation between temperature and brightness, / oc 
T^. Correlating three aim defines the bispectrum of 
e, BT'l!?-"'' = gZ7r"~^hhh. where bi,i,i, is its re- 



duced bispectrum and t/™l7!^™^ are Gaunt coefficients 



I1I2I3 

The nonlinear conversion from h-^i^i^ to bi^i^i^, 



the 



usual reduced bispectrum for temperature fluctuations 
e = 5T/T, is given by eq. g. 
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FIG. 1: Reduced bispectrum of O for = 6 as a function of 
h = ^3 (black line), compared to the analytic approximation 
from the spacial coordinate rescaling ([2| (dashed-red line). 
The other lines are the different contributions to the bispec- 
trum. In the squeezed limit, i.e. for ^3 ^ 6, the agreement is 
at the few % level. 



FIG. 2: Reduced bispectrum for li = h/10 as a function of 
h = I'i- At small I the bispectrum satisfies eq. ([2]|; deviations 
around I3 ~ 100 are due to integrated effects such as vector 
and RS. For h > 50 {I3 > 500) the long mode is inside the 
Hubble radius at recombination and the analytic approxima- 
tion fails. 



The squeezed limit: Locally, the effect of a long 
wavelength mode of C can be seen as a redefinition of the 
coordinates for the short modes [22|. At recombination, 
the long mode shifts the function e^'^ in ([5| through a 
time redefinition, e^'^ — > e^'^ + e(e^'^)', where e is the 
redefinition of the time coordinate, 77 — >■ 77 -I- e(r]) |19| . in- 
duced by the long mode. Equivalently, this can be seen as 
a long- wavelength perturbation in the free-electron den- 
sity [B], 

5, + ^ = (fe)- . (8) 

The coordinate redefinition acts also on the source S in 
eq. ([6]), generating second-order nonhnearities. For in- 
stance, it acts on FasF^F + eF + (x^dtF,^ where 
the second term comes from a rescaling of the spacial co- 
ordinates, -> a;*(l -I- C) [IS]- While A and <I> behave 
as F, £' ~ 4($ + transforms as its time derivative, 
E ^ E + (eE)' + Qx^diE. Using these properties and 
eq. ^ in eq. ([6]), one sees that the long mode simply 
induces a spacial coordinate rescaling of the source 5, 

S + (x'd^S , (9) 

while the time redefinition is re-absorbed by the long 
wavelength free-electron density The spacial coordi- 
nate redefinition above yields the bispectrum of eq. ([2]). 
In Figs. [1] [2] and [3] we compare the reduced bispectrum 



In the limit where one of the mode is super-horizon, quadratic 
terms in the velocity v or products of v times A;^ with Z > in 
F and vector and tensor perturbations in E become negligible. 
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FIG. 3; The reduced bispectrum with two I fixed, I2 = I3 
as a function of h. As expected, for h < 100 the difference 
between the computed bispectrum and the analytic approxi- 
mation is of order ~ {Ii/Ih)^- 



of 0, bi^i^i^, computed using CosmoLib2 , with the an- 
alytic approximation ([2]). For this comparison we set 
17a = and use the following cosmological parameters: 
lib = 0.15, = 0.85, h = 0.6, = 0.0, A, = 2.4 • lO'^ 
and ris — 0.97. This ensures that second-order late-time 
integrated effects do not correlate with late ISW, which 
cannot be taken into account by a coordinate redefinition 
at recombination. We have split the bispectrum pro- 
duced by the code into its different second-order contri- 
butions: SW, Doppler, Rees-Sciama, vectors and tensors. 
The bispectrum is dominated by the sum of the SW and 
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FIG. 4: Signal-to-noise density, l3^^bi-^i2i^/{Ci-^Ci2Ci.^)^^'^ , for 
I3 = 1720 and li < I2 < h- The signal is peaked on squeezed 
triangles and the acoustic oscillations are clearly displayed. 
For h > 100 the amplitude and shape agree with that com- 
puted in [15] . 



Doppler effects, which cancel themselves due to acoustic 
oscillations. Where applicable, the agreement between 
the computed bispectrum and the analytic approxima- 
tion is remarkable. Residual few per cent deviations are 
mainly imputable to the inaccuracy of the 3-d integrator 
projecting the bispectrum in harmonic space. 

Shape and amplitude of the bispectrum: We 

now numerically compute the bispectrum from recombi- 
nation and compare it to a local bispectrum [55] . We use 
WMAP7 parameters: ilt = 0.0449, = 0.222, h = 0.71, 
= 0.088, As = 2.43 • 10"^ and = 0.963. Since we 
want to concentrate on the signal generated at recom- 
bination, we do not consider second-order late-time ef- 
fects such as vector, tensor and RS, that could correlate 
with the late ISW. As an illustration, in Fig. [4] we plot 
^khh/i^h^h^is)^^^ with I3 — 1720 as a function of li 
and ^2- 

To study the observability of the signal, we consider 
a full-sky cosmic variance limited experiment between 2 
and Zmax, without accounting for lensing. The Fisher 
matrix between two bispectra i and j is given by |26| 



F, 



^] — 

2<ii<i2<'3<'„ 



Ci^Ci^Ci^Ai^i^i^ 



(10) 



with Ai^i^i^ 



[3 — 1,2,6 for triangles with no, two or 
three equal sides. We find F,.ec,rcc 0.225 (0.504), 
^locioc = 0.038 (0.088) and Fioc,ioc = 0.046 (0.069) for 
^max — 2000 (2500). The signal-to-noise ratio is given by 



S/N = F}il,, = 0.47 (0.71) ; 



(11) 



most of it is concentrated in the squeezed configuration 
so that eq. ^ and the local correction ([T]) capture most 



of the signal.^ Indeed, the contamination to a local signal 
agrees with [T9ll20] .3 



/i 



loc 
NL 



loc.i 



loc.lo 



0.82 (1.27) , 



(12) 



and will negligibly contaminate Planck's searches for lo- 
cal non-Gaussianity. 

Discussion: We have presented the most complete 
calculation of second-order effects at recombination. It 
confirms that the projection on the local shape from re- 
combination is well described by using the analytic esti- 
mate given in p!UH^ and suggests that claims of a large 
local contamination may be biased by numerical errors. 
Our calculation did not include lensing and other late- 
time integrated effects which correlate with the linear 
ISW effect. Since, for instance, lensing at recombina- 
tion partially cancels with the rescaling of the coordinates 
[THl [in] , we must include all these effects to correctly es- 
timate the total bias and the observability of the signal. 
We have shown that CosmoLib2'^'^ is a roboust tool to 
do this and, more generally, study second-order pertur- 
bations in cosmology. 
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